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FREQUENCIES OF A BRANCHED TIMOSHENKO BEAM 
by Francis J. Shaker 
Lewis Research Center 

SUMMARY 

A method is presented in this report for calculating the normal modes and frequen- 
cies of a branched Timoshenko beam. The method is essentially a modified Stodola 
method and requires an iteration procedure to determine the normal modes and frequen- 
cies of the system. In this method, an arbitrary deflection, consistent with the boundary 
condition, is assumed. Also, because of the presence of the spring-mass system, the 
frequency is not a constant factor in the governing equations and must also be assumed. 
Knowing the deflection and frequency, the corresponding shears, moments, slopes, and 
new deflection can be determined by integrating the governing differential equations of 
the system. The new deflection is then adjusted to satisfy the boundary conditions, and 
a new frequency is subsequently calculated. The process is then repeated, and the fre- 
quency is used as the criterion for convergence. Although the method can be applied for 
any type of boundary conditions, particular attention was given to those boundary condi- 
tions of interest in launch-vehicle dynamics, namely, the free-free beam and the 
cantilevered-free beam. 

In general, the iteration routine will always converge to the lowest or fundamental 
mode of vibration. To determine the higher modes, the lower mode components are re- 
moved from the assumed shear and moment distribution by utilizing the second orthogo- 
nality condition of the branched beam. The iteration routine will then converge to the 
lowest mode whose components are not removed. 


INTRODUCTION 

Experimental and theoretical studies have shown that a launch vehicle airframe can 
be idealized as a nonuniform beam for purposes of studying its lateral bending dynamics. 
When major components of the vehicle are cantilevered from the airframe, the vehicle 



can be represented as a beam with cantilevered branch beams. With such a representa- 
tion, the bending dynamics of the components and their effect on gross vehicle dynamics 
can be studied. For example, the engine shroud or the engine thrust structure could be 
idealized as a branched beam at the aft portion of the vehicle, while the payload or the 
nose fairing could be represented as a branched beam at the forward portion. For multi- 
staged vehicles, the engine and engine- support structure frequently are suspended in the 
interstage- adapter well, and this component can be represented as an intermediate 
branch. 

Because of the low aspect ratio (i. e. , L/D) of current vehicles, the bending dynamics 
of airframes are accurately described by beam theory, only if the effects of shear flexi- 
bility and rotational inertia of the beam elements are included. This beam theory, known 
as the Timoshenko beam theory, is given in reference 1. 

There are several methods for determining the normal modes and frequencies of a 
nonuniform Timoshenko beam. For example, reference 2 presents a modified Myklestad 
method which takes into account both shear flexibility and rotary inertia of the beam ele- 
ments, while reference 3 gives a matrix iteration technique which accounts for rotational 
inertia. Shear flexibility could also be handled in the second method by including the 
shear deformation in the matrix of influence coefficients. Reference 4 gives a modified 
Stodola method which also includes the effects of shear flexibility and rotational inertia of 
beam elements. The effects of flexible branches could be handled by extending either of 
these methods or by using the component mode method presented in reference 5. This 
report extends the method of reference 4 to account for the effect of the branches on the 
normal modes and frequencies of the system. 


SYMBOLS 


E 

G 

I 

*cgi 

J i 


effective shear area 

j. v J.'L 

distance from tip of x branch beam to center of gravity of i tn rigid mass 

length of i branch beam 

Young's modulus of elasticity 

modulus of elasticity in shear 

bending moment of inertia 

mass moment of inertia of i rigid mass about center of gravity 


,th 


,th 


mass moment of inertia of i rigid mass about tip of i branch beam 

f Vi 

torsion spring rate attached to i in branch beam 
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I 


M bi«i> 

M ei 


m(x, t), M(x) 
m ci 

Pi 

Qbi(^) 


q(x,t),Q(x) 

q^Qi 


Y(x,t),y(x) 

Y(x,t),y(x) 

Yo^^o 

Y' Q (t),y' Q 

a i 

h 

n 

A i 

6 

77i(^) 

0 i 


length of primary beam 

moment distribution along i 11 branch beam 

i-i. 

moment transferred to end of i Ln branch beam due to motion of spring- 
mass system 

moment distribution along primary beam 

mass of rigid body attached to 1 1 branch beam 

moment transferred to primary beam due to motion of i tn beam 

natural frequency of i 1 spring-mass system 

i-t_ 

shear distribution along i n branch beam 

ii. 

shear transferred to end of i in branch beam due to motion of spring-mass 
system 

shear distribution along primary beam 

ii. 

vertical shear transferred to primary beam due to motion of i Ln branch 
beam 

longitudinal coordinate of primary beam 

longitudinal distance from the end of primary beam to attachment point of 
i Ln branch beam 

vertical deflection function of primary beam 
bending deflection of primary beam assumed fixed at x = 0 
vertical deflection of primary beam at x = 0 
bending slope of primary beam at x = 0 

total angle of rotation of spring-mass system relative to horizontal 

1L 

bending slope at tip of i Ln branch beam 

th 

bending slope of primary beam at attachment of i branch beam 

th 

vertical deflection at tip of i branch beam 

th 

vertical deflection of primary beam at attachment station of i branch 
beam 

th 

vertical deflection function of l branch beam 

th 

bending deflection of i branch beam assumed fixed at attachment station 

angle of rotation of spring-mass system relative to static equilibrium po- 
sition 
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P A 

?(x,t),<p(x) 

$(x) 

HQ 

01 

Subscripts: 
bi, ei, i 
j,k, P, r, v 


+Vi 

generalized mass corresponding to j normal mode 

th 

longitudinal coordinate of i branch beam 
mass per unit length 
bending slope of primary beam 
shear deflection of primary beam 
bending slope of i branch beam 
circular frequency 

integers referring to branch beams 
integers referring to normal modes 


t 




THEORETICAL ANALYSIS 

Consider a launch vehicle idealized as a nonuniform beam with a cantilevered branch 
beam attached at station x . as shown in figure 1. To make the analysis more general by 
showing how a spring-mass component can be handled, assume that a rigid body mass is 
attached to the end of the branch beam by a pin connection with a torsion spring whose 
spring rate is K-. The governing equations will first be developed for a primary beam 
with one branch, and then the equations will be generalized to include (in principle) any 
number of branches. 



Figure 1. - Idealized launch vehicle coordinate system. 


4 


I 





Primary Beam Equations 


An element taken along the length of the primary beam is shown in figure 2. The 
concentrated moment and shear transferred to the primary beam from the branch beam 
are represented on the differential element by m^6(x - x.)dx and q.6(x - x.)dx, respec- 
tively. The expression 6(x - x^) is the Dirac delta function. Some of the properties of 
the Dirac delta function that will be used in the subsequent analysis are as follows: 


6(x - x.)dx = h(x - x.) = J 


X < Xj 


X > Xj 


( 1 ) 


f f 5(x - x.)dx dx = f h(x - x,)dx = «| 

Jo Jo 1 Jo 


x < X. 


( 2 ) 


x - Xj x > Xj 
v. 


where h(x - Xj) is the unit step function. The equations of motion of the primary beam 
are determined by summing forces and moments acting on the element and neglecting dif- 
ferentials of higher order. Thus, for forces, 


9q ( X> + q^x, t)6(x - x.) = -pAY(x, t) (3) 

ax 1 1 



Figure 2. - Forces and moments acting on typical primary beam element. 
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and for moments, 


8m(x>t) - q(x, t) - m.6(x - x.) = lp$(x, t) (4) 

9x 11 

Note that the rotatary inertia term in equations (1) to (4) depends only on the bending 
slope i>, since the shear forces produce only distortion of the cross section and no rota- 
tion. From the geometry of the system (fig. 1), the bending slope is given by 

¥(x,t) = 9Y (M) +Y '(t) (5) 

ax u 

where Y^(t) is the bending slope at x = 0 and dY/dx is the bending slope of the beam 
cantilevered at x = 0. The total slope will be the bending slope minus the shear slope. 
The shear slope is subtracted from the bending slope in this case because a positive 
shear force tends to decrease the total slope as shown in figure 2; therefore, 

d . y . fe’i) = t) + YJ>(t) - /3(x, t) (6) 

ax ax u 

From elementary beam theory, the angle of shear /3 at the neutral axis is given by 

/3(x,t)=aM (7) 

AG 


where A is the effective shear area. Substituting equation (6) into equation (7) gives 

s 


9Y(x,t) = 3Y(x,t) + y , (t) _ q(x, t) 
ax ax 0 ag 

fa 


( 8 ) 


The moment- curvature relations for the Bernoulli- Euler beam and the Timoshenko beam 
are the same; therefore, 


m(x, t) = El iSfa!) = El (9) 

To eliminate the time dependency from equations (2) to (5), (8), and (9), only harmonic 
motions are considered; that is, it is assumed that the shear, moment, slope, and 
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deflection functions are of the form 


m(x,t) - M(x)e lcot , q^t) = Q.e ia,t 

q(x,t) = Q(x)e 1Wt , Y(x,t) = y(x)e 1Cut ► 
m.(t) = M.e ia,t , ¥(x,t) = cp(x)e ia>t 


( 10 ) 


From equations (3) to (6) and (8) to (10), the following ordinary differential equations are 
obtained: 


d< ^ X ) = o> 2 pAy(x) - Q.6(x - x.) 
dx 1 1 


= Q(x) + M,6(x - x.) - c c?Tp<p(x) 
dx 

*>(*) = + y b 

dx 

dy(x) = dy(x) + y , _ QM 
dx dx 0 A g G 


M(x) = El - y ( x ) - 
dx 2 

J 


(11a) 


Equations (11a) represent the required differential equations of the primary beam. In the 
Stodola method, a deflection function y(x) is assumed and equations (lla) are integrated 
to obtain a new improved deflection. The constants of integration are adjusted such that 
the boundary conditions are satisfied. The process is then repeated until it converges. 
The values of M. and Q. in these equations depend on the dynamics of the branch beam. 


Boundary Conditions For Primary Beam 

The boundary conditions for two cases will be considered: (1) a free-free beam and 
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(2) a free-cantilevered beam. Other types of boundary conditions can be handled in an 
analogous manner. 

Free-free beam. - For a beam free at x = 0 and x = Z, the shears and moments at 
the boundaries must vanish. Thus, 

Q(0) = Q(Z) = 0 (12) 


M(0) = M(Z) = 0 

The first of equations (11a) and (12) yield the relation 

-Z 


(13) 


u> 


f pAy(x)dx - Q = 0 

Jo 1 


(14) 


The second of equations (11a) and (13) yield the relation 

f Q(x)dx + M, - u> 2 f ]p(p(x)dx = 0 

Jo 1 Jo 

The first of equations (Ha) yields 


(15) 


X 

Q(x) = co 2 pAy(x)dx - Q i h(x - x.) 


(16) 


Substituting equation (16) into (15) results in 

-Z sA, 


Z x Z Z 

M. + w 2 f f pAy(x)dx dx - Q. f h(x - x,)dx - w 2 / ]p(p(x)dx(x) = 0 (17) 

1 Jo Jo Jo Jo 

It is easily shown (ref. 6) that 

f pAy(x)dx dx = f pA(l - x)y(x)dx 

Jo Jo Jo 


'0 Jo Jo 

Substituting equation (18) into (19) and rearranging yield 

~Z 


(18) 


Z co 2 /* pAy(x)dx - Q. - u> 2 f Ip(p(x)dx - w 2 f pAxy(x)dx + Q,x. + M, = 0 

Jo 1 Jo Jo 


(19) 


The first term in brackets in equation (19) is simply the shear at x = Z , which is zero for 
this case. Thus, equation (19) becomes 
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CjlJ 


/ pAxy(x)dx + o> 2 f ]p<p(x)dx - Q.x. - M. = 0 

JO ill 


( 20 ) 


Equations (14) and (20) represent the equations which must be satisfied to ensure a van- 
ishing moment and shear at the boundaries of the primary beam. 

Free-cantile ver beam . - For a beam free at x = 0 and fixed at x = i, the boundary 
conditions are 


Q(0) = 0 M(0) = 0 


rd) = -§&- y(i) = o 

A s G 


(21) 


Branch Beam Equations 


The governing differential equations of the branch beam can be developed in a manner 
analogous to those of the primary beam. The coordinate system of the branch is as shown 
in figure 1, and an element along the length is as shown in figure 2, except that the con- 
centrated shears and moments q^ and m. will not be present. The governing equations 
of the branch, corresponding to equations (11a) of the primary beam, will therefore be 
given by 


d Qbi^i> 

dCi 


= a,2 (pA) bi ?7 i (^ i ) 


dM bi^i> 

d *i 


= QbA) - « 2 aw b *w 


<S?i(£ t ) 

Vt = - + 

d *i 


dr?i(^) dr? i (^ i ) Q bi (|.) 

y \ 


d ?i d« t 1 (A g G) 


bi 


M bi = (El) b . 


d\(^) 


( 22 ) 
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where the subscripts bi refer to the i^ 1 branch (for one branch i = 1) and y^ is the 
bending slope at the attachment station x.. 


Branch Beam Boundary Conditions 

To obtain expressions for the boundary conditions of the branch beam, the free-body 
diagram of the branch beam shown in figure 3 is considered. In keeping with beam theory, 
the slopes and y. are assumed to be small (i. e. , tan /3. = /3. and tan y. = y.). 

For small deflections, Newton’s second law gives the following equations for the 
spring-mass system in harmonic motion; 


or 


> 



(23) 


M ei = K i 0 i = w [ J i< 0 i + *i> + m ci a i A il 


(24) 


where J^, the mass moment of inertia of the spring-mass system about the tip of the 
branch beam, may be expressed as 



Figure 3. - Coordinate system for i^ branch beam. 
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Solving equation (24) for 0. yields 



( 25 ) 


where Pj is the natural frequency of the simple spring- mass system defined by 



Substituting equation (25) into equations (23) and (24) yields the following expressions for 
shear and moment at the right side of the branch beam. 


Q • = 

^ei 


-co 2 m ci 


1 - 


u> 


[l-^-k + a A 

K, 


r i 


M e , = KjO, = 


or 




(J./3. + m ci a. A t ) 


(26) 


(27) 


The shear and moment at the left side of the branch beam can be obtained by integrating 
the first and second of equations (22). Thus, 


<3bi (b i> = " 2 jC * < ' ,A) bl T! l <{ l >di i + «bi<°> 


(28) 


and 


M b i<»i> = Q »i (? i )d5 i - “ 2 I s flWbf l<«i> d «i + M bi<°> 


(29) 
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M i = -" 2 1/ bi <e A >bi^i(^i + / bi (*>bi^i> d *i 


L 



These two equations represent the shear and moment transferred to the primary beam as 
a result of the dynamics of the branch-beam and spring-mass system. 


Orthogonality Conditions 

To determine the orthogonality conditions for the branched Timoshenko beam, elimi- 
nate m and q from equations (3) and (4) by using equations (8) and (9) to obtain 


— (El j + A G[— - <f> ) - m.6(x - x.) = IpS 
9x V dxj s \dx ' 1 1 


d_ 

dx 


A s G (--^ 
s 'ax ) 


- q i 6(x - x t ) = pAY 


> 


(37) 


1L 

For a system vibrating in the j normal mode, 


Y(x,t) = yj( x )e 1Ct, ^ t , m.(t) = M i;j e 






$(x,t) = ^j(x)e lwt , q^t) = Q y e 


a ict> jt 


j 


(38) 
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From equations (37) and (38), 


(EIc pj)’ + [A s G(y! - ^)] - M..6(x - x.) = -w*I p( p. 
[ A s G ( y j " ' %j 6 ( x ' x i> = ~ w f pAy j 


where 


y i ‘ y i (x >’ ^ 

and the primes denote d/dx. Let y k and cp k represent the deflection and slope of the 
beam vibrating in the k** 1 normal mode. Then, multiplying the first of equations (39) by 
cp k dx and the second by y k dx, integrating over the length, and adding the resulting 
equations yields 

-a;* jT l (pA yj y k + lpcp ]( p k ) dx + M ij/ik + = jf l (El *))' cp k dx 


A s G ( y ] ■ y k ^ + [ A s G ( y ] ' ^j^k dx ] ( 40 ) 


Integrating the first two integrals on the right side of equation (40) by parts yields 


£ (f> Ay j y t + ♦ M ljnk + QySjj. = [Elvj«> k ] 0 - £ dx 


+ [ A s G(yt - ?>j)y k ]' - £ A s G(yj - ^ dx + £ A s G(y! - ^ dx (41) 


Now, for conventional boundary conditions (i. e. , free, pinned, or fixed), the terms in 
brackets in equation (41) vanish. Thus, 
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■"? X + + + Q ij 6 ik = " ^ EL ^k <* 


- X A s G(y j " ^j )(y k ' ^k^ < 42 > 

Following the same procedure for the branch beam, the equation corresponding to equa- 
tion (41) will be 

/ b ‘ t^bi-Vik + = [(“w*]* 1 * [( A .oy»Sj - f'ijhttij 1 

■ Z * 1 (E 1 W'if , 'ik d «i - / b ‘ < A S G >bl ( ”li - - ^lk> d «i < 43 > 


Evaluating the terms in brackets 


b. 




[^Vi^ikl/^e^ik-f-Mifik) 


b. 


-[(A G) W/..-T }!.)h-i 1 1= -(Q •• A., -Q..6.,) 
LV s ; bi v i] 'ij' ik J o w en lk ^ij ik ; 


(44) 


then, substituting equations (44) into equation (43) and adding the resulting equation to 
equation (42) gives 


Cl) 


<? Ay j y k + *^ k )d*+/VA) blVik - 0p> w w4 

= ^ + A S G(y j - fjXyt - V k )]dX 

*f' l^Wk ♦ < A s G ) bl (’-ij -VfrW - *ik>] d «i - »W* AQ eij 4 ik < 45 > 


Also, from equations (23) and (24) 
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M ../3., - Q .. A., = -co. 0., (J.a.. + m .a. A.-) 
erfik ^eij lk j ik v x i] ci 1 ly 


+ + m ci (A ij + a i a ijH A ik + «!«&)] ( 46 > 


From equation (24), 

K 0 

J l“ij + m ci a i a ij = “T 1 (47) 

"J 

Substituting equation (47) into equation (46) and substituting the resulting equation into 
equation (45) gives 

"?{/* ^ A ¥k + * J^ 1 f<P A >biVik + »>blVikJ d *i * I cgi“ij“ik 

+ m cl (i,j + + a i a lk>} = \ f f Ehp 'f’k * A s G(y i - y "j)( y k - «’k>] dx 

+ K iVik + _/ b ‘ * < A s G >bi<”ij * ^’("ik - ^lk)l d «i} < 48 > 

All the terms which appear in equation (48) are symmetric with respect to j and k ex- 

O i-L 

cept cof. Therefore, if the system is considered to be vibrating in the k n mode and the 

same procedure is followed as for the j** 1 mode, an equation identical to equation (48) is 
2 2 

obtained except cof' will be replaced by co£; that is, 


co. 


[{jf °“ Ay k y J + +j[ bi t^biVik + + 'cgi^ij^ik 

+ m ci (4y + ajayXA^ + A [EI. ?J ' S ^ + A,.G(y; - ^Hy,' - <P k )]ifa 

+ K iVik + / b ‘ K*Wlk + < A S G >bi<”ij - ^ik ~ ' f 'ik^ d «i| < 49 > 
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2 2 

Subtracting equation (49) from equation (48) and assuming co. j= oof gives 

J K 

+ + j( bi [ ^ >A) bi 1J ij T, ik + 

+ 1 cgi“«“lk + m ci (a lj + + a i“ik> " 0 < 50 > 

Equation (50) represents the first orthogonality condition of the branched beam. To obtain 

9 2 

the second orthogonality condition, multiply equation (48) and (49) by and oof, re- 
spectively, and subtract. Then, for co? * 

£ l V&v'f't + A s G b’j - VjUyfc - 

+ jC ‘ f (EI) bi^ij^ik + <VVii - - ^tk^i = 0 < 51 > 

An alternate form for the second orthogonality condition is obtained by using the relations 

M = EL Q = A s G(c p - y») 


M 


bi - %i = <Wi - ’ll) 


in equation (51). Thus, 


r^ax + r 

Jo EI Jo A = G Jo (Er) “ 


d «i 


f 1 ^bij^bik 

i (A s G >bi 1 


+ eik 1U = 0 (52) 


K i 


,2 , , ,2 


Generalized mass. - Equations (50) to (52) are valid for oof i= oof. To obtain the gen- 

J k 

eral expression for these orthogonality conditions, let 
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“i ' jC ( pAy ? + H*)* 1 + [ <pA> w”« + ““VyK 


+ ‘cgi^ij + m ci (A ij + a i a ij) 2 


(53a) 

where jUj is referred to as the generalized mass of the system vibrating in the j th nor- 
mal mode. Then, in terms of the generalized mass, the first and second orthogonality 
conditions become 




+ •cgi'Vik + m ei (a y + Vij>< a ik + VW = "j 5 


jk 


(54a) 


and 


£ [ (EI ''’S«i + A s G(y i - - ''"k>l dx + K i d ij d sk 


£' [< EI >biVik + <W>y - ^ip^ik - ^ik>3 d «l - "jfyjk < 5 5a> 


where 



is the Kronecker delta function defined by 


fl 



j = k 


j *k 


Equations (54a) and (55a) satisfy equation (48) for j = k and equations (50) and (51) for 
j * k. 

Natural frequen cy. - An expression for the natural frequency of the system can be 
obtained by setting j = k in equation (48). Thus, 
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«? = 


/?>■/' 


M 


bij + Q bij 


K EI )bi ( A s G )bi 


d^ i + 


M 


ei] 


K. 


jC ( pAy i + *‘*l) dK + jT* + (Wbi^itK + I cgi“« + m ci (A Ij + a l“«> 2 


(56a) 


Primary Beam with Several Branches 


The equations developed on the preceding pages are for a primary beam with one 
branch. The results can be immediately generalized to include any number of branches 

by replacing Q.6(x - Xj) and M^x - x.) with ^ Q i 6(x - x^) and ^M.6(x - x.), respec- 

i i 

tively, in the equilibrium equations. The development will then proceed in the same man- 
ner, and the resulting equations will be similar to those previously developed, except that 
the branch beam parameters will be summed over all the branches. For example, as- 
sume that there are n branches. Then, a summary of the pertinent equations will be as 
follows. 


Governing differential equations of primary beam: 

= co 2 pAy(x) - Y Q.5(x - x.) 
dx Z-j 1 1 

i=l 

n 

MJS = Q( x ) + V M,6(x - x.) - c o 2 lpcp(x) 

h 

dx 

dy(x) _ dy(x) + y , _ Q(x) 
dx dx 0 A g G 

M(x) = EI d2 r& 

j. 2 


(lib) 



The last three equations are independent of the number of branches and are included 
herein only for later reference. 

Free-free boundary conditions: 


w 


r l 

1 J pAy(x)dx - > , Q i = 0 
i=l 


l rl n n 

CD 2 f pAxy(x)dx + u> 2 / Ip<p(x)dx - Y] Q.x. - Y* M. = 0 

y ° y ° 


>- (53b) 


Generalize mass: 


<*, ■/ K ♦ * E /* K.4 ♦ 

i=l 


E [w“« + m ci (A ii - v/l (53b) 


i=l 


Orthogonality conditions : 


l n A>. 

(pAyjy k + 1p<Pf k )<& + X 1 t^bi^ifik + ^bi^ik^i 


tl 


+ [I cgi a ij a ik + m ci (A ij + a i a ij )(A ik + a i A ik>] = ^j 6 jk < 54b) 
i=l 


n 


f [EIc p'f'^ + A g G(y! - ^(y^. - «^ k )]dx + ^ Vij 0 ik 
0 i=l 


JL x-b 


£ 1 [(EDbi^ij^ik + ( A s G ) bi ^ij - ^ij>^ik - *ik>] d ^i = w j ^j 6 ik ( 55b > 

i=l 
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Frequency equation: 




/ 


5 + i ldx+ 


EI A g Gj 


£{/ 


M b« + 


P^bi < A s G > bi 


df, + 


< 


K. 


/ (pAy* + S>,>*)dx + ^j/* 1 [lpA) bl 7)(j + Gp)^ y] d 5i + I cgl “ ij + m c i( 4 i + a i“ i]> 

i=l L 


(56b) 


MODIFIED STODOLA METHOD 

In the previous section, the governing differential equations of a branched 
Timoshenko beam have been developed. Also, the boundary conditions for the free-free 
and the free-cantilevered primary beam have been expressed in integral equation form. 

In this section, a modified Stodola method will be applied to solve the differential equa- 
tions of motion subject to the free-free and free-cantilevered type boundary conditions. 

The Stodola method, as applied to beam vibration, starts by assuming a deflected 
shape of the beam which satisfies the boundary condition. Based on this assumed deflec- 
tion, the corresponding shears, moments, slopes, and new deflections can be determined 
by integrating the governing differential equations. Then, by using the boundary condition 
equations, the new deflection is adjusted to satisfy the boundary conditions. At this point, 
the natural frequency of vibration can be determined. Using the new improved deflection, 
the process is repeated until the natural frequency converges. If the branched beam has 
an attached spring mass, as in the case considered in this report, it is also necessary to 
assume a frequency of vibration at the start of each iteration. The reason for this is that 
the frequency does not appear as a constant factor in the governing equations as developed. 
This point will be clarified as the Stodola method is applied in detail. 

For the free-free case, the iteration routine is started by assuming that the deflected 
shape of the primary beam is a polynomial of the form 


y(x) = Yp 0) ( x ) 


V=D 



(57) 


I 
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where the superscript (o) refers to the first iteration and the subscript p referes to the 
mode being determined. For a cantilevered beam, clamped at x = l, the initial assumed 
primary beam deflection can be taken as 


P 

y(x) = Y^(x) = (x - l) 

v=l 

A deflection must also be assumed for each branch beam. For the 
fleeted shape can be taken as 


x - 


vl 


(58) 

J-U 

i Ln branch, the de- 


4?^ = 6 i o) + r i o) ^i + 


where 6^ is the deflection of the primary beam at the attachment station and 


,(°) = Y'(°)(x. 

l p v l 


(59) 


is the slope at the attachment station. 

Before starting the iteration routine, it will be convenient to normalize the deflection 

to unity at the tip of the primary beam. The branch beams should be multiplied by the 

2 

same normalizing factor as the primary beam. Also for convenience, co , which appears 
as a factor, is set equal to one, since the mode shapes can only be determined to within a 
constant factor. The to terms which do not appear as a factor are set equal to one in 
the first iteration and then for subsequent iterations, the value determined by the previous 
iteration is used. 

j.r 

First, consider the iteration routine of the i Ln branch beam. For the assumed shape, 
the shear and moment at the tip of the branch beam is given by equations (26) and (27) as 


-m 


Q • = 

^ei 


ci 


i-s! 


2 

Pi 





(60) 


M e,= 


i.«d 


(J.B. + m .a. A.) 
'ri Cl 1 v 


(61) 
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Then, the shear and moment at the attachment station will be given by equations (33) and 
(34), respectively. 


Q * = ' jC ‘ + «ei 


(62) 


M i = -f / 151 (pAlbi^iW, + j£ l W, 


+ Q .b. - M 
mi i e 


(63) 


Since or /3. is not known in the first iteration, it is set equal to zero, and for subse- 
quent iterations the values from the previous iteration are used. 

The shear and moment distribution along the branch beam can now be calculated by 
integrating the first and second of equations (22) and noting that, for the adopted sign con- 
vention, 


QbifO) = Qj, M bl (0) = -M, 


(64) 


Therefore, 


Qbi(li) = f h <P A^SjJdJj + Qi 


( 65 ) 


- X 1 Q bi ( «i )d| i ■ X h - M t 


( 66 ) 


When the shear and moment distribution is known, a new improved deflection of the 
branch beam can be determined by utilizing the third, fourth, and fifth of equations (22); 
thus, 


v[(^) = 


/' 


i M bi (^) 


(El) 


d ^i + r i~ 


bi 


gbigj) 

(A s G >bi 


(67) 


and 
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To determine the values of 6. and y. in equations (67) and (68), the deflection and slope 
of the primary beam at the attachment point of the i branch must be determined. This 
can be done by an integration procedure similar to that used for the branched beam. 

For the primary beam, the shear and moment distribution can be determined by in- 
tegrating the first and second of equations (lib). Thus, 


Q(x) 


/-x 

= J pAy(x)dx - J2 Qi h ( x " x i) 
i=l 


(69) 


and 


n 


n 


M(x) = XX pAy(x)dx dx - J' Tpcp(x)<ix - h(x - x^)dx + E M t h(x - x.) 

(70) 


i=l 


i=l 


where h(x - x.) and f h(x - x,)dx are given by equations (1) and (2), respectively. 
1 Jo 1 

Then, from the third, fourth, and fifth of equations (lib), 


y(x) = 




Q(x) 

A s G 


dx + y^x + y Q 


(72) 


The constants of integration yy and y Q must now be adjusted to satisfy the boundary 
conditions. First, for compactness, the following quantities are defined: 


Y(x) = 



dx 


dx 


$(x) = 



dx 
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M bifej> 

<ED bl 


d «i d «i 



Using these substitutions, the deflection and bending slope of the beam can be written as 

y(x) = y(x) - $(x) + yJjX + y Q 


cp(x) = y’(x) + yjj 

V^i) = ~ ?i(^) + [y f (x.) + yytj + [y(x.) - $(x.) + y|jX. + y Q ] 


y 


^i(?i) =T?-(? i ) + [y'Cxj) +y{j] 


(73) 


j 


For the free-free beam, the boundary constants yfc and y 0 must be adjusted to satisfy 
equations (14) and (20). The expressions for and Mj, which should be used in these 
equations, are given by equations (35). Therefore, the following equations must be sat- 
isfied 


/ 

Jo 


pAy(x)dx + 







pAxy(x)dx + 



Jpcp(x)6x + tf (pA) bi (x. + + tf 

i=l i=l 


y ( 74 ) 


+ 





( x i + b i) 


A i + 


n 



i=l 
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If equations (73) are substituted into equations (74), two equations in two unknown con- 
stants are obtained, namely, yfo and y^. Note that in equations (74) the expressions 
for A i and / 3 . are given by 


= ^fri) = Vfti) ~ ^(bj) + [y’(x.) + y^Jbj + y( Xj ) - $(x.) + yJjX. + y Q 

^i = ^jO^) = vfai) + y ’(xj) + yj) 


( 75 ) 


Substituting equations (73) and (75) into (74) yields 


n 


f pA[y(x) - 4>(x) + yjjx + y Q ]dx 
*/0 


HI .a. 

h-fl^) + y’(xj) + yjj] 


i=l 


tJL 1 <^ A >bi{^i^i) - ?i(*i) +[y’( x i) +yb]^i + y( x i) - *( x i) + yb x i +y o} d? i 

i=l 



( 76 ) 
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PM y( x ) - ^(x) + y^x + y 0 ]dx lp[y(x) + y^dx + ^ (lp) bi [r 7[(q) + y'( x i ) + yfjd «- + ^ ^^pAfx. + - C^q) 



[y’(x t ) +y[j]l i + y(x i ) - *(x.) + y^ + y^d^ + \ — l - ( x t + b P + a i |{r7 i (b i ) - qO^) + [y'(x i ) + y[j]b. + y(x.) - $(x t ) +y|)X i +y 0 ] 


«? V j,p? 



- t m ci a i( x i + b i ^ + + y'( x i) + yol = 0 


Collecting terms in and and rearranging yields equations (78a) and (79a). 
Shear boundary condition: 


jT 1 pA[y(x) - <fr(x)]dx + ^ j^ 1 


<P A ) bi E^i^i) - CiftiWi + 



« i=i 


I'^i) + X y,(x) 1 X 1 ^bih d *i + ^ + a i> " 



i . 

,2 \ K. 


- ^(bj) + y(x t ) - $(x.)] + ^ [y< x i> “ 1 ^ A) bi d? i] 


AjC pAxdx+ i / bi <p A) M (x 1 + « i> d « i + 


^-Xi + b i + ai .^ l ( b 1+} 


m . / I„„,w \ 

— ei_ i -_£E_ = o 

,.,2 \ K , 
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Moment boundary condition: 


£ pM y(x)-$(x)]dx + ^ Jpy’(x)dx + Z jC l (pA)bi<Xl + + E jC' 



[y( x i) - *( x i)] JT 1 (P A W x i + *i )d *i + j — f m ci a i (x i + b i + a i> + 



+ 


AjC 


pAx dx + 


E /** ^ A) bi (x i + + 

i=l 



x i + a i +b i-^( x i +b i )u ’ 2 
K i 


J 


(79a) 
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To put equations (78a) and (79a) into a more convenient form, suitable for calcula- 
tions, we define the constants given in the appendix. The subscripts r which appear 
with the constants K in this appendix signify that K depends on the particular mode be- 
ing sought. The other constants are independent of the mode and depend only on the mass 
and geometric properties of the system. Using the appendix, equations (78a) and (79a) 
can be written as 


A lr y 0 + A 2r y 0 = C lr 

(78b) 

A 3r y 0 + A lr y 0 = C 2r 

(79b) 


where 
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n 

C lr = K lr + E (K lri ' K 4ri + K 5ri + Vl ' K 10ri> 
i=l 

C 2r = K 2r ‘ K 3r + £ (K 2ri ' K 3ri ' K 6ri + ^ri + K Sri ' K llri> 
i=l 

Solving equations (78b) and (79b) yields, for yfo and y^, 


C lr 

A 2r 


A lr 

C lr 

C 2r 

A lr 

V = - 

A 3r 

C 2r 

A lr 

A 2r 

y 0 

A lr 

A„ 

2r 

A 3r 

A lr 


A 3r 

A i 

lr 


(80) 


Equations (80) give the boundary constants for the free-free case. For the free- 
cantilevered case, we obtain from equations (21) the following expressions for yfo 
and y Q : 


yq = -y'(i) y 0 = $(D - y (D + 1 • y '(i) (8i) 

After the integration constants have been determined from either equations (80) or (81), 
a new deflection can then be obtained from equations (73). The frequency of vibration can 
now be calculated from equation (56b). Before starting the next iteration, it will be nec- 
essary, in general, to recalculate yQ and yjj for the free-free case based on the new 
improved frequency. This is required because the assumed deflection must satisfy the 
boundary conditions. Then, after determining a new improved deflection, the deflection 
is normalized, and the process is repeated until two successive frequencies are equal. 

A flow chart of the iteration routine is given in figure 4. 
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Free-cantilevered 










































High Modes of Vibration 


The iteration routine described in the previous section will always converge to the 
first or fundamental mode of vibration. To determine the higher modes of vibration, it is 
necessary to remove all lower mode components from the assumed shape. Let y*(x) and 
V* i(£) represent the assumed shape of the r Ln mode. Assume that the mode shapes of 
all the v = r - 1 modes are known and that it is required to determine the r m mode. 
Since any arbitrary shape can be represented by a series involving the normal modes of 
the system. 


y*(x) = Cjy^x) + C 2 y 2 (x) + . . . + C v y v (x) + . . ." 

r?*i(£) = + c 2 T ?2i^ ) + • • • + c v %i ( £) + • • • 

J 


(82) 


where the C y 's are unknown constants that must be determined. The shear and moment 
distribution, corresponding to the assumed deflection, can be expressed by the series 


Q*(x) = CjQjCx) + C 2 Q 2 (x) + . . . + C y Q v (x) + . 
M*(x) = CjMjCx) + C 2 M 2 (x) + . . . + C v M y (x) + . 
Q ri (0 = CjQ^) + C 2 Q 2 .(^) + . . . + C v Q vi (£) + . 
M r .(4) = + C 2 M 2 .a) + . . . + C v M yi (|) + 




>- 


(83) 


For | = bj, the last equation becomes 


M* . 
eri 


“ C l M eli + C 


2 M e2i + 


CyM 


evi 


(84) 


Next, multiply the first of equations (83) by [Q v (x)/A q G]dx, the second by 
[M y (x)/EI]dx, the third by [Q vi (^)/( A s G ) bL d ^ and the fourth by [M vi (4)/(EI) bi ]d£., and 
integrate over their respective lengths. Then multiply equation (84) by Mg^/PL. 
These operations yield the following set of equations: 
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/ 




dx 


AG 


= C, 


f l QiQ v ** 

i 


/ 

</o 


M*M dx 
r v 

El 


= C, 


f 1 M l M v 
/ El 

*/o 


dx 


+ • 


f lQ?iQ y_ i d^ = c 1 f 1 QliQvi 

7 ^ 7 (AsG) m 


d|+ . 


/' 

«/o 


M ri M vi 

( E D b i 


dl = Cj 


/‘ 

*/o 


( EI )bi 


d£ + 



(85) 


M* .M . M -.M • 

err evi _ eli evi , 

— L/ + • 


K i 


K i 


C v M evi 

K i 


+ • 


Adding equations (85) and using the second orthogonality condition which is given by equa- 
tion (55b) yield 



M*M Q*Q 
- — - + — — -|dx + 
El AG 

o 



Q*.Q . 
^rrvi 


( A s G ) 


bi 




M* .M . 
en evi 

K 4 


v^v 


(86) 


where is the natural frequency of the v Ln mode and /j . is the generalized mass of 

th v v 

the v mode given by equation (53b). 

For n branches, equation (86) becomes 
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( 87 ) 



M r M v , Qr Q v 
El A g G 


dx + 



M*.M . Q*.Q . 
rx VI. TTX^Vl 

T 


(ED 


bi 


(A S G). 


bi 


M* -M . 
d , + _erx_evxl 

K i 


0 ) 2 JU 

v p v 


Knowing a y , the purified shear and moment distribution becomes 


q-1 

= «? - E C v Q v 
V=1 

q-1 

M r = M r ~ X) a v M v 
v=l 

q-1 

q . = q*. . y c q . 

^rx rrx / y v wi 
v=l 

q-1 

M ri = M*. - V C.M., 


^ 1 v* vi 

V=1 


( 88 ) 


Equations (88) represent the shear and moment distribution which must be used in the 
iteration routine to determine the higher modes of vibration (i. e. , 2, 3, . . . )• 


CONCLUDING REMARKS 

A computer program, based on the analysis presented in this report, has been writ- 
ten for finding the normal modes and frequencies of a branched Timoshenko beam. Sev- 
eral test cases have been run and the results obtained show that the iteration routine con- 
verges in a reasonable length of time. A flow chart showing the iteration routine is given 
in figure 4. An examination of this flow chart shows that a number of integrals must be 
evaluated during each iteration. Any of the standard methods of numerical integration 
can be used to evaluate these integrals, and the best method will depend, in general, on 
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the mass and stiffness distribution along the beam. The summation method, which is 
thoroughly covered in reference 7, is particularly suitable for most practical cases. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 23, 1967, 

180-06-06-03-22. 
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APPENDIX - DEFINITION OF CONSTANTS USED IN EQUATIONS (80) 



N, 


= f pAx dx 

*/0 


n m ■ f b ' <p A )w< x i + 


N ci = m cl< x i + b i + a i> 


M, 


= f l pAx 2 


dx 


M bi = <P A >bi< x i + d? i 


M ci = m ci< x i + b i + “i>‘ 



P bi = / «”bi^ 
M di = m ci (x i + b i )2 
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K lr = j[ l 

K lrl - <pA> W l?l(?l) ' ’’ l(5l)ld * 1 

K 2r = Jq 1 pM^( x ) • yC*)]* 

K 2ri = J[ 1 ^ A W x i + 

S'l __ 

K 3r = / lpy'(x)dx 

**ri - jC l VhfflW 
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K 6 ri = y , ( X i)1 M bi- N bi X i +P bi + 


( \ 

1 


1 - 


CO 


V A) 


I M 

M . - N -X. + I . - -Hii — (M,. - N,.x.) 
ci ci i cgi v di di r 


Kjr i = W - y( x i)] N bi 




K 8ri ^ 


^ N ci - N di iw - + - y< x i)] 


i - 


CO 2 \ ~ K i 


K 9ri = ™ bi [$(x) " yW] 


K 10ri = 


/ \ 


m ci a i 


1 - — 


v pf/ 




K llri " 




l-d 


Ph 
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